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A probabilistic structural analysis of a space shuttle main engine (SSME) 
turbopump blade has been underway at the Lewis Research Center for the last 
3 years. The immediate objectives of this study are to evaluate the effects 
of actual variations, also called uncertainties, in geometry and material prop- 
erties on the structural response of the turbopump blades. In this study a 
normal distribution has been assumed to represent the uncertainties statisti- 
cally. Uncertainties were assumed to be totally random, partially correlated, 
and fully correlated. The magnitudes of these uncertainties have been repre- 
sented in terms of mean and variance. These two quantities were selected such 
that the absolute magnitudes of uncertainty at any point were either less than 
or equal to 10 percent of their original values. 

Many studies have demonstrated that probabilistic analysis methods for compo- 
nents under random loading are more reliable than deterministic approaches. 
Probabilistic methods have been predominantly used in fatigue, fracture mechan- 
ics, and structural reliability analyses under ramdon vibrations (e.g., stud- 
ies reported by Dover, 1979; Yang, 1981; Wirshing, 1981; Kawamoto, 1982; and 
Huang and Nagpal , 1984). Particularly in fatigue, improvements in estimating 
fatigue life ranged up to several hundred percent in comparison with a widely 
used deterministic approach developed by Miner (1945). Because of the poten- 
tial, the application of probabilistic methods has been accepted in the other 
fields of engineering in recent years. These fields include input loading, fi- 
nite elements, and metallurgy. A probabilistic approach to develop a compos- 
ite loading for SSME components is underway at the Battelle Laboratories under 
the technical guidance of R. Kurth (1985) and has been discussed by Shinozuka 
and Lin (1981). Other studies using probabilistic finite-element methods have 
been briefly reviewed by DasGupta (1986). A variational approach for develop- 
ing the probabilistic finite elements is under development by Belytschko and 
Liu (1985). The usefulness and importance of the probabilistic approach, espe- 
cially for turbopump blades, has been summarized by Chamis (1986). 

Blade response, recorded in terms of displacements, natural frequencies, and 
maximum stress, was evaluated and plotted in the form of probabilistic distri- 
butions under combined uncertainties. These distributions provide an estimate 
of the range of magnitudes of the response and probability of occurrence of a 
given response. Most importantly, these distributions provide the information 
needed to estimate quantitatively the risk in a structural design. 

*Work performed for the Structural Mechanics Branch under contract NAS3. 

2-219 


PRECEDING PAGE BLANK NOT FILMED 


A wide range of response, such as a maximum stress distribution, implies that 
smaLl natural uncertainties in the geometry and material properties can cause a 
large variation in the response. Consequently, designs based on deterministic 
methods will not be as safe as probability-based designs. Higher stresses, 
which have low but certain probability of occurrence, are not included in the 
deterministic designs. In most cases available experimental results are ^ 
insufficient to provide an entire range of stresses and their probabilities of 
occurrence. Limited experimental results as well as deterministic methods pro- 
vide estimates of the mean response but do not cover the scatter. Consequent- 
ly, the risk in design even after using the factor of safety may or may not be 
within acceptable limits. 

The results of the study indicate that an uncertainty up to 10 percent in the 
material properties had no significant effect on the response. However, uncer- 
tainty up to 10 percent in geometry, only through the thickness, has signifi- 
cant effect on the structural response. Based on the results, a number of 
probabilistic models has also been developed using regression analysis. These 
models predict the response for given uncertainties in the geometry and mate- 
rial property parameters. The models for uncertainties in temperature com- 
bined with those in geometry and material properties are under development. 



OBJECTIVE 


The analysis was performed on the airfoil of a typical SSME turbine blade. 
Nominal dimensions of the airfoil are 1 to I /2 inches long, about 1 inch wide, 
and a & inch thick. Generally, these blades are made of anisotropic materi- 
als. The number of blades per turbine stage is about 65. The other parts of 
the blade are the platform, shank, and fir tree. 
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PROBABILISTIC STRUCTURAL RESPONSE 

The probabilistic structural response of the blade was recorded in terms of 
three natural frequencies, root stress, which is also the maximum 
stress, and blade tip displacements. The root stress is the stress at the 
junction of the airfoil and the platform. Since the blade airfoil is consid- 
ered fixed at this junction for this analysis, stress at that location becomes 
the maximum. The tip displacements were estimated in three directions. Two 
displacements were in the plane of the airfoil, and one was perpendicular to 
the plane of the airfoil. 
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methodology 


The developed methodology is generic and can be used to analyze other 
structural components which may or may not be related to this problem. This 
chart describes the major steps used for performing this analysis. The first 
step is the generation of random numbers. Correlated or uncorrelated random 
numbers are used for perturbation of the geometry, spatial location, and/or ma- 
terial properties. Uncorrelated numbers indicate no bearing of the perturba- 
tion of one point on any other point on the airfoil. 

The second step is modeling the blade with a finite-element model. Eighty tri- 
angular shell elements with 55 nodes were used to model the airfoil. 

The third step involves setting up of the design to run simulations. Full fac- 
torial design was used because it is considered more economical than other 
available approaches. 

The fourth step involved analyzing the response and developing its probabilis- 
tic distribution and models to predict the response for given uncertainties. 
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SAMPLE RANDOMNESS USED IN PERTURBATIONS 


Both uncorrelated and correlated random numbers are generated by the comput- 
er. They can be generated using any distribution and with preselected statis- 
tical properties such as mean, standard deviation, etc. In the present study 
a normal distribution with selected mean and variance was used. The random 
numbers have the property that they can be multiplied by any constant without 
losing their randomness. A plot of random numbers is shown. 
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FRACTORIAL DESIGN FOR THREE VARIABLES 


A factorial design is an economical method for performing real experiments or 
computational simulations. Unlike a parametric approach, the influence of 
individual variables and their interactions can be estimated discretely. 

Another advantage of the factorial design is that the range of the variables 
can be extended by amending the design scientifically. 

There are two limitations with the factorial design: it can be used only for 

the independent variables, and the results of the study are only applicable 
within the range studied. 

A representation of the factorial design for three variables is shown in the 
attached figure. The solid circles at the corners or in the center of the box 
represent the experiment or simulation points. A range of variables is defined 
as the difference between its two extreme levels. 



VARIABLE 1 
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RANDOM COORDINATE PERTURBATION OF SSME BLADE 


A geometric perturbation is done by perturbing the nodal coordinates in all 
three dimensions. The figure shows an exaggeration of perturbation of nodal 
coordinates in the x direction only. Perturbation in all three directions 
is difficult to show in graphical form. The geometric perturbation in this 
study is intended to simulate realistic geometric variation due to manufactur- 
ing tolerances and the errors in finite-element modeling. The variations in 
geometry at any point were limited to less than 10 percent. 

Perturbations were created using random numbers. Uncorrelated and partially 
correlated random numbers were used to simulate truly uncorrelated or patterned 
variations in the geometry. 
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RANDOM YOUNG'S MODULUS DISTRIBUTIONS IN SSME BLADE 


Material properties perturbation is done by perturbing the properties of each 
finite element. The attached figure is a demonstration of the variation of mo- 
dulus of elasticity for a blade that is made of an isotopic material. For a 
blade made of anisotropic material the entire material property matrix is per- 
turbed at the same time. Material properties, like geometry, were perturbed 
using uncorrelated and partially correlated random numbers. 

These perturbations are to represent realistic variation in the material prop- 
erties. The maximum variation at any location was limited to 10 percent. The 
"10 percent" number was selected based on expert opinion. 
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PROBABILISTIC MODELS 


Probabilistic models were developed using regression analysis. The general 
form of the model is represented at the top of the following figure. The left 
side of the model represents dependent variables, which can be any variable in 
the response. The right side of the model consists of means and standard devi- 
ations of uncertainties and their coefficients. Over 20 models were developed. 
For demonstration purposes coefficients of only one of the models are plotted. 
This model is for estimating the variation first natural frequency due to com- 
bined effect of uncertainties in geometric and material properties. 


MODEL 

VMMBLE 17 Y = ( C0NST ) + (C0EF) 1 <r, + . . . + {C0EF) 6 a C33 + . . . + (C0EF) 21 a 3 a Cs3 
WHERE 

ffj STANDARD DEVIATION OF PERTURBATION IN I DIRECTION 

oCjj STANDARD DEVIATION OF PERTURBATIONS OF MATERIAL PROPERTY MATRIX 
ELEMENT Cjj 
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INFLUENCE OF MATERIAL PERTURBATIONS ON NATURAL-FREQUENCY DISTRIBUTION 


The response variables were evaluated for both uncorrelated and correlated ran- 
dom variations in geometry and material properties. The variation in the natu- 
ral frequencies for partially correlated variations in material properties was 
slightly wider than for uncorrelated variations. However, the increment in 
variations was not significant. The response for variation in geometry alone 
and combined variation in geometry and material properties was studied. 

Sample histograms of variations in the first three natural frequencies are 
shown . 


UNCORRELATED PARTIALLY CORRELATED 



THIRD 


NATURAL FREQUENCIES, kHz CD-88-32640 


2-230 


• METHODOLOGY TO QUANTIFY UNCERTAINTIES IN THE DESCRIPTORS OF SSME BLADES 
HAS BEEN DEVELOPED AND APPLIED. 

• METHODOLOGY IS GENERIC AND EXTENDABLE TO ANY STRUCTURAL COMPONENT AND 
ALL ASPECTS OF STRUCTURAL ANALYSIS AND DESIGN. 

• GEOMETRIC UNCERTAINTIES SHOWED SIGNIFICANT EFFECT ON COMPONENT 
STRUCTURAL RESPONSE. 

• MATERIAL PROPERTIES UNCERTAINTIES AND THEIR INTERACTIONS HAVE 
INSIGNIFICANT EFFECTS ON COMPONENT STRUCTURAL RESPONSE. 

• RANGE OF VARIATION IN STRUCTURAL RESPONSE QUANTIFIED AND EXPLICIT 
PROBABILISTIC MODELS HAVE BEEN DEVELOPED. 

• INFLUENCE OF UNCORRELATED AND PARTIALLY CORRELATED VARIATIONS IN 
GEOMETRY AND MATERIAL PROPERTIES ON THE RESPONSE HAVE BEEN STUDIED. 
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